Phononic self-energy effects and superconductivity in CaC6 



A. Sanna 1 ' 2 , S. Pittalis 3 , J. K. Dewhurst 1 - 2 , M. Monni 4 , 
S. Sharma 1 ' 2 , G. Ummarino 5 , S. Massidda 4 and E. K. U. Gross 1,2 
1 M ax- Planck- Institut fur Mikrostrukturphysik, Weinberg 2, D-06120 Halle (Germany) 
2 European Theoretical Spectroscopy Facility (ETSF) 
3 Department of Physics and Astronomy, University of Missouri-Columbia, Columbia MO, 65211 (USA) 
4 Dipartimento di Scienze Fisiche, Universita degli Studi di Cagliari, Cagliari (Italy) and 
C.N.I.S.M. Dipartimento di Fisica Politecnico di Torino, 
C.so Duca degli Abruzzi 24, 10129 Torino (Italy) 

We study the graphite intercalated compound CaC6 by means of Eliashberg theory. We perform 
an analysis of the electron-phonon coupling and define a minimal 6-band anisotropic structure, that 
leads to a Fermi surface dependance of the superconducting gap. A comparison of the superconduct- 
ing gap structure obtained using the Eliashberg and the superconducting density functional theory 
is performed. We further report the anisotropic properties of the electronic spectral function, the 
polaronic quasi-particle branches and their interplay with Bogoljubov excitations. 

PACS numbers: 74.25.Jb, 74.25.Kc,74.70.Ad 



5h . 

C/3 



The electron-phonon interaction leads to many signif- 
icant physical phenomenon in solids (notably, supercon- 
ductivity), and has therefore been studied extensively 
both in model systems and in real materials. One impor- 
tant aspect of this kind of interaction is the formation of 
a coupled electron-phonon system with new interesting 
features such as the appearance of polaronic sub-bands 
branching from the main electronic bands. This low en- 
ergy features of the electronic structure can be observed 
thanks to the recent developments in the resolution of 
Angle- Resolved Photo Emission Spectroscopy (ARPES). 

The theoretical background to deal with metallic po- 
larons has been laid down by Engclsberg and Schricf- 
fer (ES)i. For this they used a field theoretical ap- 
proach combined with Einstein-Debye model to mimic 
the phonon spectrum. The ES theory shows the damping 
of electrons by phonons and the development of branches 
in the electronic dispersion corresponding to energy and 
strength of the phonon modes. In superconductors, the 
electron-phonon interaction leads to the formation of 
a superconducting gap below the critical temperature 
T c . These phenomena can be well described within the 
Eliashberg theory which extends the ES theory to the 
superconducting stated and reduces to the ES theory in 
the non-superconducting normal state. 

Due to several computational complexities, a proper 
account of the material specific electronic and phononic 
structures could not be achieved until very recently- 
Eiguren et al£ and Eiguren and Ambrosch-Draxl 3 - stud- 
ied the effect of the electron-phonon interaction on the 
electronic self-energy in the normal state. The main 
properties of the spectral function in the superconduct- 
ing state have been reported by Scalapino^ and, more 
recently have also been studied using the ARPES exper- 
imental data£~— . However, to the best of our knowledge, 
no first-principles attempt has been made to study the 
effect of polarons in the superconducting state. 

In the present work we use the Eliashber g 4 i 5 i 23 i 24 
method to study the behavior of ES polarons; a detailed 



analysis of the electronic self-energy, including electron- 
phonon contributions is performed. In particular, the 
features originating from the anisotropy of the electron- 
phonon coupling are investigated. Most importantly, it 
is shown how the polaronic branches change in the su- 
perconducting state below T c . 

The system that is considered for this analysis is the 
graphite intercalated compound CaCg This material has 
the highest superconducting T c observed so far (11.5 K) 
amongst the group of graphitic compounds. Graphite 
related materialas have attracted considerable interest in 
the last few years, mostly due to the appealing possibil- 
ity of tuning their physical properties^. In particular 
it is possible to vary the conductivity of graphite from 
semi-metallioii to metallic and to superconducting^— 
by adjusting the level of intercalation. In Ca interca- 
lated graphite superconductivity arises from the strong 
electron-phonon coupling provided both by C and Ca 
phonon mode a 16 ' 21 . This coupling is strongly anisotropic 
with C and Ca related phonons acting selectively on 
the multiple Fermi Surface (FS) sheets of the system 2 ^. 
These peculiarities make the system particularly inter- 
esting. 

The paper is organized as follows: In Sec. I, the main 
concepts and physical quantities describing our results 
are introduced by reviewing the Eliashberg theory of su- 
perconductivity. Sec. II is devoted to a detailed descrip- 
tion of the principal computational techniques employed 
in this work. In Sec. Ill, results for CaCg are discussed. 
Sec. IIIA reports on the structure of the electron-phonon 
interaction. Sec. IIIB focuses on the numerical solutions 
of the Eliashberg equations with a //* determined using 
results from Density Functional Theory for Superconduc- 
tors (SCDFT) calculations. In Sec. IIIC, the polaronic 
features of the excitation spectrum of CaCg are eluci- 
dated, both in the normal state (Sec. IIIC1) and in the 
superconducting state (Sec. IIIC2). Finally, conclusions 
arc drawn in Sec. IV. 
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I. METHODS 

Central quantity in the Nambu-Gor'kov formalism of 
superconductivity is the 2x2 Green's functional: 



G(k, iu! r , 



G(k,iu n ) -F(k,iuj n ) 
-F*(k,icj n ) G(k,-iu n ) 



(1) 



where G(k, itu n ) and F(k, iu> n ) are, respectively, the nor- 
mal and anomalous electronic Green's functions in re- 
ciprocal space. u n are the Fermionic Matsubara fre- 
quencies given by ui n = ?r(2n + l)k B T with T be- 
ing the temperature and k B the Boltzmann constant. 
Following a well established procedure^ 3 -, non-interating 
Kohn-Sham system with Green's function Go(k, iuj n ) — 
[ibj n UQ — CkO"3] , is used as a starting point. Here 
<Tj (j = ... 3) are the Pauli matrices and £k are the 
Kohn-Sham eigenvalues relative to the Fermi energy. The 
interacting Green's function can then be obtained using 
perturbation theory: 



[G(k,icj n )] 1 = [G (k, 



IUJ„ 



£(k,iu; n ). (2) 



The following approximation for the electronic self- 
energy is used: 



Efctw,,) = -k B T J2 ^ 3 G(k', 
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Matsubara integration to a low energy ( a fraction of cV 
) by a renormalization procedure ^ — ► fi* introduced by 
Morel and Anderson^ ' 27 ' 31 , fj,* can be calculated within 
the random phase approximatio n 22 ' 28 ' 29 (RPA) and the 
resulting T c are usually in reasonable agreement with ex- 
periments^ 3 -^. However, usually //* is adjusted to obtain 
the experimental T c . 

The self-energy in Eq. ^ is k-dependent and leads to 
anisotropic Eliashberg equations which are computation- 
ally very demanding. A simplification can be introduced 
retaining a minimal anisotropic structure needed for the 
properties of interest. The FS may be divided into por- 
tions (FS sheets), with each sheet identifying a corre- 
sponding intersecting energy band. These FS sheets and 
energy bands can be labeled using the same index (say, 
J). We shall refer to such a division as to a multi-band 
decomposition (details of our procedure are given in Sec. 

MM- 

The electron-phonon coupling can be averaged over a 
prescribed multi-band divisions (see below, Eq.s (fTOj) and 
(jlip ). Corresponding, multi-band resolved self-energy, 
Sj, expanded in the basis of Pauli matrices has a form: 

£j(rt) = iui n [1 - Zj{n)] a + Aj(n)Zj(n)ai , (4) 

(terms which only result in a rigid shift of the Fermi level 
are neglected). The multi-band resolved Green's function 
reads: 



^ |ffk,kv| 2 A/(k - k',iuj„ - iu n ') 
+ criW(k-k')] (3) 

here D is the phonon propagator (D = 
— 2w qiI //(o; 2 + w 2 jly )), 3k,k',y are the electron-phonon 
matrix elements 2 ^ between states with wavevector k and 
k' , and due to a phonon mode of index v and wavevector 
q=k — k' and w q .„ is the frequency of the mode obtained 
via linear response 2 ^ of the Kohn-Sham system. This 
way of calculating w q! „ is known to lead to a very good 
agreement with the measured phononic branches, at 
least for standard metals and insulators 2 ^. W in Eq. 
[3] is the screened static electron-electron interaction, 
it accounts for those parts of the interaction which do 
not involve any phononic contribution. The a\ factor 
in front of W accounts for the fact that exchange and 
correlation effects arc already included in Go(k, iu> n )2&. 
Then only off-diagonal contributions of W are retained 
in S. 

The treatment of the Coulomb term needs particular 
care. Within Eliashberg theory, an arbitrary cut-off in 
k space is needed in order to avoid serious convergency 
problems in the Matsubara summational^ 7 -. A conven- 
tional way to deal with this problem is to choose an en- 
ergy cut-off of the order of the Fermi energy, and to as- 
sume that the product of W with the density of states 
(DOS) is constant: fj, = WN, N being the DOS per spin 
at the Fermi energy. It is then possible to restrict the 



Gj(k, iu n ) = 



iw n Zj{n) + £j k Aj(n)Zj(n) 
Aj(n)Zj(n) iuj n Zj (n) - Ok 

[Zj(n)u; n } 2 + + [A,,{n)Zj{n)} 2 



(5) 

where £,/k are the Kohn-Sham eigenvalues in the J-th 
band. Using Eqs. and (|5|) in the Dyson equation, 
we arrive at the following set of coupled self-consistent 
equations^i^ 3 .! 3 ^: 

Zj{n) = 1 + — > \jj,{n,m) ,(6) 

«w n ~, Rj'(m) 



A p j (n) = irk B T ^ \jj,(n,m) 



771, J' 



Aj>(m)Zji(m) 
Rj>(m) 



,(7) 



\C ( \ i ^sr * Aj>(m)Zj> (to) 
a j W = -Kk B T Hj r 5-7^ , (8) 



m, J' 



Rj>(m) 



Rj(n) = ^/(co 2 +A 2 (n))Z 2 (n) , 
2Loa 2 Fjj,(Lu) 



Xjji (n, to) = / du) 



(9) 
(10) 



Here Aj{n) = Af{n)+A C j{n) is the total superconduct- 
ing gap accounting for phononic and Coulombic contri- 
butions on the J-th FS sheet at frequency icj n . Zj(n) is 
the (phononic) mass renormalization function. This term 
enters in the diagonal part of the electronic self-energy 
and contributes both to the superconducting state and 
the normal state. Due to the assumption that all the 
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diagonal contributions stemming from W are already ac- 
counted at the level of the normal state Kohn-Sham sys- 
tem, Zj(n) has a purely phononic character. 

The off-diagonal Coulombic contributions are ac- 
counted by the FS-dependent . In this work, we 
choose in such a way to reproduce the gap struc- 
ture obtained within SCDF T 22 ' 29 (we shall come back to 
this point in Sec. lIII B|) . Since A ph (n) goes to zero for fre- 
quencies much larger than the phononic scale, the cut off 
u! c in the Matsubara frequency introduced for Coulom- 
bic terms can be uniformly applied to all the terms of the 
Eliashberg equations. 

a 2 Fjj/(uj) in Eq. (fTUj) is the band-resolved Eliashberg 
function. This quantity results from the averaging of the 
electron-phonon interaction over FS sheets: 



a 2 Fjj,(u) 



1 



E 

keJ 
k'ej' 



(11) 

where Nj is the density of states per spin at the J-th FS 
portion. 

Once the gap function Aj(n) and the mass renormal- 
ization function Zj(n) have been obtained on the imag- 
inary axis solving the Eliashberg equations (Eqs. 
(HUl)), they can be efficiently continued to the real axis 
via Pade approximant tcchnique22r— , allowing the com- 
putation of the real-axis retarded Green's function. In 
particular, we will deal with G := G 11 which is given by: 



ujZj(uj) + £ 



7k 



[uZj{u)Y -ej±-[&j{u)Zj{u)Y 



(12) 



with the corresponding spectral function defined as: 



1 



ImGj(k,w). 



(13) 



The energy and life time of a quasi-particle are given 
by the real part and (the absolute value of the) imagi- 
nary part of the pole positions, respectively. The spectral 
function may have broader structures and, thus, the con- 
cept of quasi-particle may not apply rigorously. In order 
to extract the main features, we look for the zero, z p , of 
the denominator of G, i.e. we look for solutions of the 
equation: 



Zj(z p ) 



(14) 



in the complex plane. 



II. COMPUTATIONAL DETAILS 

Electronic eigenvalues, phononinc frequencies, and 
electron-phonon matrix elements are calculated using the 
ESPRESSO pscudopotential based package^ 2 .. All cal- 
culations are done using the GGA (Generalized- Gradient 



Approximation) with the Perdcw-WangH parameteriza- 
tion for the exchange-correlation functional. Ultrasoft 
pscudopotentials 4 ^ are employed. A 30 Ry cut-off is fixed 
for the planewave expansion of the wavefunctions and 300 
Ry for the electronic charge. The Brillouin zone is sam- 
pled with a 6 x 6 x 6 k-point grid, and electron-phonon 
matrix elements are obtained on a 10 x 10 x 10 grid. More 
details can be found in Ref. [22|. 

The double Brillouin zone integration appearing in 
the definition of the band-resolved Eliashberg functions 
in Eq. (ITT1) is evaluated with a Metropolis integration 
scheme. Random k-points are generated on the Brillouin 
zone, then each k-point is accepted or rejected with a 
probability depending on £ k , and its weight is set in- 
versely proportional to the acceptance probability. We 
use a set of about 2 ■ 10 4 accepted k-points per band. 
Then electron-phonon matrix elements on this random 
mesh are obtained via interpolation from those calculated 
on the regular grid. 

Eliashberg equations are solved using a 2500 mcV cut- 
off of the Matsubara frequencies, with the cut-off ui c for 
the Coulombic interaction equal to 500 meV. Both pa- 
rameters are much larger than the maximum phonon fre- 
quency of CaCg, which is about 200 meV. The number 
of Matsubara frequencies at each temperature is fixed 
by the energy cut-offs, and the M Matsubara frequen- 
cies on the positive imaginary axis are used to construct 
(M, M) Pade approximants which are used for the ana- 
lytic continuations^ to the real axis. 



III. RESULTS AND DISCUSSION 

A. Properties of the electron-phonon coupling 

We calculate the Eliashberg function defined by 
Eq. (|11[) in three different ways: 

• Averaging the electron-phonon coupling over the 
full surface: referred to as 1-band FS or isotropic 
approximation. 

• Splitting the FS into three parts: this is shown with 
three different colors in the FS (see Table U). The di- 
vision of the Fermi surface leads to a division in the 
electronic bands shown in the band structure plot 
of Fig. [5j an electronic band and a portion of the 
FS have same color if they intersect. The first por- 
tion FS 1 is the external FS sheet (shown in green), 
which comes from 7r states. FS 2 (shown in blue) 
is the spherical Ca Fermi surface. FS 3 (shown in 
red) is the 7r-prism, a two-dimensional FS having 
the shape of a hexagonal prism which crosses the 
spherical Ca Fermi surface (the corresponding elec- 
tronic band is the band 3). This division is referred 
to as the 3-band FS approximation. 

• Splitting the FS into six parts: each of the three 
portions in the 3-band FS approximation is further 
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Energy (meV) 

FIG. 1: (color online) 3-band FS resolved a 2 F,,j,(uj) of 
CaC6. The three panels refer to the incoming FS, from which 
the electron is scattered, while the three lines in each panel 
refer to the outgoing FS. 



split into two parts. The above ir (external) FS has 
been divided into a less coupled^ outer part with 
\k xy \ > 0.4 a.u., and the rest; these two portions 
are referred to as la and lb portion, respectively. 
The Ca spherical Fermi surface is cut into 2a por- 
tion (with |k z | < 0.18 a.u.) and the 2b portion 
(the rest). 3a and 3b portions for the w- prism are 
defined in a similar manner. The same boundaries 
are used to further split the corresponding energy 
bands. This overall division is called the 6-band FS 
approximation. 

In Tabic U is presented the DOS at the Fermi level, the 
intra- and the inter-band electron-phonon couplings: 

Xjji = 2 / a2FjAuj) du> . (15) 
J w 

In the 3-band FS approximation the interaction is 
dominated by the off-diagonal coupling terms, especially 
by the inter-band scattering from states on the spherical 
Ca Fermi surface to states on the 7r-bands. The main rea- 
son for this is a strong electron-phonon coupling in the 
former (from 2 to 1) and large DOS in the latter (from 
2 to 3). The full 3-band a 2 Fjj,(uj) matrix (see Fig. [J) 
shows the distribution of the coupling among the various 
phonon modes- band 2 couples strongly with Ca modes 
(giving a low frequency peak around 10 meV) while, band 
1 couples mainly with the high-frequency stretching C 
modes (at 170 meV). Band 3 shows the most homo- 
geneous coupling (its intra-band spectral function looks 
similar in shape to the total Eliashbcrg function). 

The further decomposition into the 6-band FS approx- 
imation doesn't introduce qualitative differences with 
respect to the 3-band decomposition but, as we shall 
show, it results in a better quantitative description of 
the anisotropy of the superconducting properties. 
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TABLE I: Electron-phonon coupling Xjji and DOS for the 
3-band FS (top) and 6-band FS (bottom) divisions. DOS is 
given in states/eV/spin.The first index runs over the columns. 
Isotropic DOS and A are 0.787 and 0.870, respectively. The 
picture on the top represents the FS of CaCe: the external 
portion in green; the internal sphere in blue; and the hexagon 
that cuts the sphere in red. The labels "a" and "b" correspond 
to the additional splitting defining the 6-band division. 



B. Solution of Eliashberg equations 

The Eliashberg equations are solved in three different 
ways corresponding to three ways in which the Eliash- 
berg function defined in section Mil The solution to the 
Eliashberg equations lead to a strong anisotropy in the 
gap (see Fig. [5]) ; the smallest gap corresponds to the ex- 
ternal FS (FS 1), while the highest value of the gap is 
related to the 2a structure which forms the central part 
of the Ca spherical Fermi surface. We note that at the 
phononic level, the anisotropic structure obtained here 
agrees very well with the one obtained within SCDFT 
in Ref. if (see panels (c) and (f) in Fig. EJ. If the T c 
is determined from this gap function, without including 
the Coulomb interaction, it is not strongly affected by the 
multi-band character (i.e. the anisotropy of the gap func- 
tion); the isotropic T c is about 33.5 K and only slightly 
higher in the 6-band case with a value of 34 K. 

In order to include the Coulomb interactions the 
matrix fijj, is needed. Typically is determined 

by fitting to the experimental data. However, for 
CaC6 the SCDFT gap well reproduces the experimen- 
tal measurements^— . Therefore, in this work, we de- 
termine /i* 7 j, fitting to the gap structure obtained from 
SCDFT. 

Interestingly, the simple semi-isotropic approximation: 
H*j j, = /i*Nj/N, turns out to be sufficient to reproduce 
the SCDFT gap. In the 6-band FS approximation the ex- 
perimental T c ~ 11.5 K is reproduced for n* ~ 0.21 [see 
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FIG. 2: (color online) (a) Anisotropy of the superconducting 
gap, in the 3-band FS approximation; and (b) 6-band FS ap- 
proximation. In (a), the dark green thick line is used for the 
gap computed over FS 1; the blue thin line for the gap com- 
puted over FS 2; and the red dashed line for gap computed over 
FS 3. In (b), dark green thick and dotted thin lines are used 
for la and lb FS, respectively; blue thin and dash-dotted thick 
lines for 2a and 2b FS, respectively; and red dashed thick and 
dashed thin lines for 3a and 3b FS, respectively, (a) and (b) 
refer to phonon only calculations Inset (c) shows the effect of 
the inclusion of the Coulombic interaction: the Eliashberg T c is 
given as a function of fi* (see text for details). Inset (d) reports 
the low-temperature multi-gap structure for fi* — 0.21. Panels 
(e) and (f) report the gap distribution functions (phonon only 
calculations) in the 6-band FS approximation in SCDFT— and 
Eliashberg theory (with a gaussian broadening applied), respec- 
tively. 
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FIG. 3: (color online) Eliashberg functions, (a) Supercon- 
ducting gap and mass renormalization function calculated 
solving the isotropic Eliashberg equations on the imaginary 
axis (b) Analytically continued gap function on the real axis, 
(c) Analytically continued mass renormalization function on 
the real axis. Full lines are used for the real part and dashed 
lines for the imaginary part of these complex functions. On 
the imaginary axis the functions are purely real. 



interaction, A c (n) is both independent of k and the Mat- 
subara frequencies. More physical features emerge from 
the analytic continuation to the real axis. In particular, 
we sec a three-peak structure in both the real and imag- 
inary parts of Z(n) that correlates with the peaks in the 
a 2 Fj. r (Lu). 



C. Analysis of self-energy effects 



inset (c) in Fig. [2]. With this choice of /i} 7 , the inclusion 
of Coulombic effects reduces the T c without significantly 
affecting the anisotropic structure of the superconducting 
gap. The only difference is that by including the Coulom- 
bic interaction the gap corresponding to the 2b portion 
of the FS becomes slightly larger than the 3a portion 
(see Fig. [2Id)). This choice of a semi-isotropic Coulom- 
bic pairing which is often done in Eliashberg theory can 
be validated by the present analysis. However, this can- 
not be applied as a general rule; it has been shown, in 
cases like in MgB ^ 46 ' 47 , that a more detailed Coulombic 
structure is necessary to get the correct gap anisotropy. 

Both A(n) and Z(n) are purely real- valued on the Mat- 
subara frequencies. For the isotropic case the frequency 
dependence of A(n) and Z(n) is shown in Fig. [SJa). It 
is clear that Z{n) has a value of 1 + A for small \n\ and 
then monotonically decreases to 1 at energies much larger 
than the available phonons. This behavior is almost inde- 
pendent of the values of fi* and temperature. The A(n) 
function also monotonically decreases as a function of in- 
creasing energy. The low-energy value is the fundamen- 
tal superconducting gap, while in the high-energy limit 
approaches A c (n). For an isotropic and static Coulomb 



1. Normal state: isotropic approximation 

We discuss in this section the simplest case in which 
the Eliashberg equations are solved above T c . Since T c is 
quite low with respect to the phonon energies, solution 
above T c is almost equivalent to imposing the solution 
at A = when T = 0. In this section, we use only 
the isotropic solutions of the Eliashberg equations and a 
parabolic band dispersion. 

The solution of the Eliashberg equations give Z(ui) 
which is shown in Fig. (3{c) (black line). Using Z(lj) as 
input, Eq. (fT4")) is solved in order to obtain the quasi- 
particle dispersion curves reported in Fig. |4ja). This 
figure shows how the unperturbed electronic band, by 
getting dressed with the phononic self energy, develop 
branches in correspondence with the three main phononic 
peaks in the Eliashberg function. These electronic quasi- 
particles dressed by the electron-phonon interaction, so 
called polarons, are nearly dispersionless. Only one dis- 
persive branch, that goes to zero from about 50 meV, 
is observed. This structure also has a very short life- 
time (~ 100 meV). All the other polaronic modes have 
instead a longer lifetime (between 1 and 10 meV) and 
thus appear as sharp quasi-particles. However, the po- 
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laronic branches, carry very little of the total spectral 
weight. Most of the spectral weight is still localized near 
the bare electron dispersion line. This can also be seen 
in Fig. HJb) , where the spectral function is shown in the 
same energy/momentum window as the quasi-particle 
plot. In this case, we see how the main electron band 
acquires a finite lifetime and instead of branchings only 
kinks appear. These kinks correspond to the three main 
peaks in the a 2 F(ui). 

To move from a qualitative description to a more quan- 
titative one, the spectral function A(w) is examined (see 
Fig. Etc)). At k = kp more than 50% of spectral weight 
is accounted for by a single peak of infinite lifetime at the 
Fermi energy. This peak (single green line in Fig. HJc)) 
moves to higher energies with increasing k- vector (slowly 
growing in width) up to an energy of about 10 meV, 
where it merges with the polaronic branch generated 
by the low- frequency Ca modes (blue long-dashed line). 
Above 10 meV the peak is broader (blue thick line) be- 
cause the electrons can relax through the generation of 
Ca phonons. This broad peak then behaves in a similar 
way as the narrow peak below 10 meV; i.e., it increases in 
energy with k up until it merges with another polaronic 
band which originates from the low-frequency C modes 
(red dot-dashed line) and has an energy ~ 50 meV. It 
becomes very broad (yellow short-dashed line) and is dif- 
ficult to follow as it merges with the high-frequency C 
mode. This behavior is very similar to the case of the 
Einstein phonons discussed by Engelsberg and Schrief- 
fer in Ref. |l|. In our case, this is due to the three-peak 
structure of the a 2 F(uj) i.e. due to the combined effect of 
the two-dimensionality of graphite along with the pres- 
ence of weakly bound Ca ions. Therefore, at the isotropic 
level the self-energy effects in CaCg (type described by 
Eq. ([3]) ) are particularly simple. 



2. Normal and superconducting states: anisotropic features 

The degree of complexity of our analysis is increased by 
making use of the real KS band dispersions of CaC e 16 ' 22 . 
Multiple FSs that couple with different phonon branches 
are accounted for by adopting the 6-band decomposition. 
As shown in Fig. [SJc), (d), (e), and (f), the Ca band 
couples mostly with low-frequency modes, therefore it 
shows polaronic structures only up to 50 meV. One single 
kink can be observed in the spectral function at 50 meV. 
The kink at 10 meV is not visible simply because below 
this frequency the spectral function itself is just a sharp 
peak. 

The band that has most structures is the one that pro- 
duces the 7r- prism FS (Fig. Efg) , (h), and (i)), due to 
the coupling of all three sets of modes. The external 
FS, which couples mostly with high-frequency C modes, 
shows only a weak kink around 160 meV. The pola- 
ronic branchings (Fig. 03V), (c'), (d'), (h'), and (i')) 
have similar structures as observed in the isotropic limit. 
Anisotropic features are less marked than in the spec- 



tral function, because Eq. (TT4)) used to determine these 
features does not retain information about the spectral 
weight in the branches. 

The spectral function in the superconducting phase is 
gapped and the multi-gap features can be clearly seen in 
the lower panels of Fig. Mb*), (c*), (P), and (g*). The 
gap ranges from 1.5 meV around the point b in the band 
structure (Fig. [5jb*)) to about 2.2 meV around point c 
(Fig. Etc*)). The spectral function shows a (textbook- 
like) hyperbolic dispersion- this is the signature of Bo- 
goljubov excitations 2 ^, as compared to the normal elec- 
tronic excitations which have a linear band dispersion 
(see Fig- EJc* ) and (c#)). The distance between the ver- 
tices of the hyperbola is equal to 2A(0). As the temper- 
ature rises to approach T c , A(0) — ^ and the hyperbola 
tends to its asymptotes. At the same time the spec- 
tral weight of the two reflected components (right part 
of the upper branch and left part of the lower branch in 
Fig-EJc*)) also goes to zero and the excitation spectrum 
becomes normal (Fig. [5jc#)). 

One of the two arms of this hyperbolic dispersion cor- 
responds to the normal electronic dispersion line, and it 
behaves in a similar way as in the non superconducting 
state. The other arm that is a unique feature of the Bo- 
goljubov excitations loses spectral weight as the distance 
from the FS increases. However as it reaches the energy 
of the Ca in-plane modes (from 10 to 15 meV) it deviates 
from the hyperbolic arm and follows a polaronic (dispcr- 
sionless) behavior. One can appreciate the formation of 
this dispersionlcss Bogoljubov polaron in the upper-right 
and lower- left corners of panel (b*) in Fig. [5] 



IV. CONCLUSIONS 

Superconductivity in the graphite intercalated com- 
pound CaC6 is studied using Eliashbcrg theory and su- 
perconducting density functional theory. Within a multi- 
band description and assuming a structureless Coulomb 
interaction, we performed a detailed analysis of the 
influence of strongly anisotropic electron-phonon cou- 
pling on the k-dependence of the superconducting gap. 
Anisotropics computed with Eliashbcrg theory and su- 
perconducting density functional theory were found to 
be in very good agreement with each other and with ex- 
periments^ . 

In this context, from the solution of the Eliashbcrg 
equations, we have shown how anisotropic polaronic 
bands emerge over different Fermi surface sheets. The 
interplay between superconducting (Bogoljubov) excita- 
tions and polarons has also been studied. We reported, 
for the first time, how Engelsberg-Schrieffer polarons 
evolve from the normal state to the superconducting state 
in CaCg. 
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FIG. 4: (color online) Polaronic bands and spectral function for the normal state in the isotropic approximation. 
On panel (a), the polaronic quasi-particle dispersion curves are shown. This is obtained from the real and imaginary 
part of the solution of Eq. (|14[) . The inverse lifetime of the state is given by \Im(z p )\ and it is reported for some 
points as "error-bar". Panel (b) shows the spectral function A(k,uj). The abscissa represent the |k| axis, and the 
point in which A(k,u) crosses is kf. In ordinate, the frequency axis is reported. The color-scale goes from zero to 
about 15 eV" 1 (note that, we cut it off to enhance the structures of A(k, ui)). Panel (c) shows A(k, ui) for different 
values of the momentum. A few values are highlighted and correspond to k- values indicated in panel (b). Panel 
(d) shows a 2 F(uj) computed in the the isotropic 1-band approximation. 
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FIG. 5: (color online) Polaronic bands and spectral function for the anisotropic and superconducting states of CaCe 
determined within the 6-band approximation. Central panel shows the band structure of CaCe: the continuous 
green line is the Cn band corresponding to FS 1; the blue dot-dashed line is the Ca inter-layer band corresponding 
to FS 2; and the red dashed line is the second Ctt band corresponding to FS 3 (see Table U for plots of the several 
FS portions). Further subdivisions are not reported explicitly. Panels from (a) to (i) show the spectral function 
near the Fermi energy in a wide energy window of 300 meV. The letters correspond to the labeling of the crossing 
points in the central panel. The color-scale goes from zero (black) to 15 eV _1 (white). Higher peaks in the spectral 
function have been cut off. Panels with a primed index show the quasi-particle spectrum (see text for details). 
These also show the polaronic branches together with the main electronic band. Colors are proportional to the 
line-width of the quasi-particle state: black corresponds to zero and yellow/white correspond to about 0.1 eV. 
Panels (b*), (c*), (f), and (h*) - focusing on the superconducting gap (a logarithmic color-scale is used) - show 
a zoom of the spectral function near the Fermi energy in a narrow energy window of 20 meV. Panel (c*) is the 
same as (c*) but evaluated in the non-superconducting state. 



